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Abstract 

We study the existence of stationary solutions of the Vlasov- 
Poisson system with finite radius and finite mass in the stellar dy- 
namics case. So far, the existence of such solutions is only known 
under the assumption of spherical symmetry. Using the implicit func- 
tion theorem we show that certain stationary, spherically symmetric 
solutions can be embedded in one parameter families of stationary, ax- 
ially symmetric solutions with finite radius and finite mass. In general, 
these new steady states have non-vanishing average velocity field, but 
they can also be constructed such that their velocity field does vanish, 
in which case they are called static. 



1 Introduction 

Large stellar systems such as galaxies or globular clusters are often described 
by a density function / = f(t,x,v ) > on phase space; t G IR denotes time and 
x,v GiR 3 denote position and velocity respectively. Under the assumption 
that the mass points in the ensemble, i. e., the stars, interact only by the 
gravitational potential which they create collectively and that in particular 
collisions are negligible, the time evolution of the ensemble is described by 
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the following nonlinear system of partial differential equations, known as the 
Vlasov-Poisson system: 

dtf+vV x f-V x U-V v f = 0, 

AU = 4irp, p(t,x) = J f(t,x,v)dv. 

Here U = U(t,x) denotes the gravitational potential of the ensemble and 
p = p(t,x) denotes its spatial mass density. For simplicity we assume that 
all particles have the same mass, equal to unity, and set the gravitational 
constant equal to unity as well. 

In the present investigation we are interested in the existence and proper- 
ties of solutions of this system, which are independent of time. Such solutions 
are usually called stationary. If they have the additional property that their 
average velocity field vanishes, i. e., Jvf(x,v)dv/p(x) = 0, xGlR 3 , we shall 
call them static. If U is independent of time, the energy 

E(x,v):=^v 2 + U(x) (1.1) 

of a particle with coordinates (x,t>)GlR 6 is constant along solutions of the 
characteristic system 

x — v, v — —VU{x) 
of the Vlasov equation. Therefore, the ansatz 

f(x,v) = <!>(E) 

automatically satisfies the Vlasov equation and reduces the problem of find- 
ing a stationary solution of the Vlasov-Poisson system to solving the semi- 
linear elliptic problem 

AU = 4nh*(U) (1.2) 

where h<& is obtained by inserting the ansatz for / into the definition of p. If 
other invariants of the characteristic flow are known — such as the modulus 
of the angular momentum \x x v\ in case of spherical symmetry — then $ can 
depend on these additional invariants as well and the right hand side of (|1.2Q 
can become explicitly x-dependent. The main difficulty with this approach 
is to show that a solution of (|1.2j) — once its existence is established — leads 
to a stationary model with physically reasonable properties, in particular, 
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with finite mass and finite radius, i. e., p is compactly supported. In P| 
this program was carried out under the assumption of spherical symmetry, 
where it can also be shown that the distribution function / must be of the 
form f(x,v) = $>(E, \x x v\) for some $. Spherically symmetric stationary 
solutions are automatically static. To the best of our knowledge, the existence 
of static or even stationary solutions of the Vlasov-Poisson system in the 
stellar dynamics case with finite radius and finite mass and without spherical 
symmetry is open. This is interesting because for a self gravitating fluid it is 
known that every static solution must be spherically symmetric, cf. [13[ . 



In the present paper we show that this is not so for a selfgravitating en- 
semble as described by the Vlasov-Poisson system. In fact, we will show that 
every static solution (f ,p ,Uo) in a certain subclass of the spherically sym- 
metric ones obtained in || is embedded in continuous one-parameter families 
(/ 7 ,p 7 ,[/ 7 ) of stationary solutions with axial symmetry which coincide with 
{fo,Po,U ) for 7 = 0, have finite radius and finite mass, and are not spheri- 
cally symmetric for 7 7^ 0. Families of static as well as families of stationary 
but not static such solutions are obtained for the same spherically symmet- 
ric steady state. For the precise statement of our result we refer to the next 
section. The basic idea of the proof is the following: Assuming that U is 
axially symmetric, i. e., U(Ax) = U (x) for every x G H 3 and every rotation A 
around, say, the X3-axis, the quantity 

P(x,v) :=XiV2-x 2 v 1 , (1.3) 

that is the ^-component of the angular momentum of the particle with co- 
ordinates (x,v) GlR 6 , is conserved along characteristics. We make the ansatz 

f(x,v) = <S>(E,jP) 

with $ such that 7 = leads to one of the spherically symmetric solutions 
with finite radius and finite mass obtained in H. Then we transform the 
problem ( |1.2| ) to the problem of finding zeros of an operator T(7, •) where 
for 7 = we know a zero, namely Uq, and we can try to prove our result by 
applying the implicit function theorem. The central idea which makes this 
approach work is to look for U 1 as a deformation of U , i. e., £7 7 (x) = U (g(x)) 
for some diffeomorphism g on IR 3 , and to formulate the problem of finding 
zeros of T over the space of such deformations instead of the space of poten- 
tials. Whereas the original problem (|1.2|) had to be solved on IR 3 , it turns out 
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that one needs to know the deformation only on a compact neighbourhood 
of the support of the original solution (fo,po,Uo), and this provides useful 
compactness properties. In particular, this deformation approach is essential 
in proving that the derivative of T at Uo is an isomorphism. Finite radius 
and finite mass of the resulting stationary solutions are then immediate con- 
sequences of the corresponding properties of {fo,Po,U ). 

The approach which we explained above has been used by Lichtenstein 
for proving the existence of slowly rotating Newtonian stars, as described by 
selfgravitating fluid balls, cf. |H| A translation of Lichtenstein 's ap- 
proach into modern mathematical language and the framework of the implicit 
function theorem is due to Heilig, cf . || , and the present paper owes much 
to that investigation. Our approach is analogous to || but the actual proofs 
are different, so that we decided to give a self-contained presentation of the 
arguments for the present case of the Vlasov-Poisson system. Our paper pro- 
ceeds as follows: In the next section we formulate our result and the general 
framework for its proof. In particular, we define the Banach spaces which 
serve as domain and range for the operator T(j,-), introduce the deforma- 
tion mappings and show how our result is obtained from the implicit function 
theorem. The continuous Frechet-differentiability of T with respect to the 
second argument and the fact that at zero this derivative is an isomorphism 
are then established in Sections 3 and 4 respectively. 

We conclude this introduction with some references to the now quite ex- 
tensive literature on the Vlasov-Poisson system. Global existence of classical 
solutions has been established in |TJ| , cf. also [10, 13, 20, 21] and the review 



article [jig] . For the plasma physics case, where the sign of the source term 
in the Poisson equation is reversed, the existence of stationary solutions, say, 
on bounded domains or with a fixed ion background or external force field, is 



much easier to obtain, cf. for example |Io |. Moreover, there are now several 
results on the stability properties of such stationary plasmas, cf. for example 
[^, |T7|. The stability question for the stellar dynamics case is much harder, 
and preliminary results can be found in |4], |22| , cf. also |19 . Coming back to 



the topic of the present paper we mention that in families of stationary 
solutions of the stellar dynamic Vlasov-Poisson system with axial symmetry 
were obtained, but these models have infinite mass and infinite radius. 



4 



2 The main result 



In this section we give the precise formulation of our result and show how it is 
obtained from the implicit function theorem, postponing the rather technical 
verification of the assumptions of the latter to the last two sections. We hope 
that most of our notation and terminology is self-explaining, but the following 
needs to be introduced: The closed ball in IR 3 with center and radius R > 
is denoted by 

B R :={xEJR 3 \ |x|<R}, 

and 

B R :=B R \{0}; 
\x\ denotes the Euclidean norm of xGlR 3 . Also, let 

Si:={:reIR 3 | |x| = 1}, 

denote the unit sphere in IR 3 , and denote the line segment joining two points 
x,x'GlR 3 by 

i^:={Ax + (l-A)x'| AG [0,1]}. 
The set of transformations which are to leave our solutions invariant is 

is a rotation around the x 3 -axis 
or the mapping IR 3 3 (xi,x 2 ,x 3 ) i-> (xi,x 2 , — x 3 )| ; 

note that in addition to axial symmetry we require reflection symmetry with 
respect to the plane {x3 = 0}. Let 

C s (B R ):={feC(B R )\f(Ax)=f(x), AeS, xeB R }. 

Clearly, 

V/(0) = 0, fzC\B R )nC a (B R ), 

and this is the reason for introducing the extra reflection symmetry. 

For the phase space distribution function / of our stationary solution we 
make the ansatz 

f(x,v):=(f>(E)if)(rfP), x,veM 3 

where 7 G IR and U is assumed to be axially symmetric; the quantities E and 
P were defined in ( |1 . 1|) and (|1.3| ). Throughout our investigation, 0:IR— ► 
[0,oo[ and ip : IR— > [0,oo[ will satisfy the following assumptions: 
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(01) G Lf oc (IR) for some p > 2, and there exists a constant E G IR such that 
0(E) = for £? > P a. e., and <p(E) > for £ < E a. e. 

(02) The ansatz f Q (x,v) =<f)(E) leads to a nontrivial, static solution 
(/oiAb^o) °f the Vlasov-Poisson system, which is spherically symmet- 
ric, i. e., po and U depend only on \x\, and such that p G C^(1R 3 ) with 
suppp = -Bi and ?7oGC 2 (IR 3 ) with Mm^^ooUo^x) = 0. 

(V>) ^eC 2 (IR) with ^'(0) = and ^(P) = 1^P = 0. 

Remark: For E G IR and — | < p < | the function 



0(E) := 



(E -EY,E<E , 
,E>E 



obviously satisfies (01) and leads to a spherically symmetric steady state 
(fo,Po,U ) with finite radius and finite mass, cf. 0, Thm. 5.4]. The solution 
has the required regularity, without loss of generality we can assume that 
supppo = -Bi, and since lim\ x \- i . 00 Uo(x) exists we can take this limit to be zero 
by redefining Eq accordingly. Thus a large class of the so-called polytropic 
steady states satisfies the assumptions (01) and (02). 

Theorem: There exists a constant 70 >0 such that for every 76]- 7o,7o[ 
there exists a nontrivial stationary solution (/ 7 ,p 7 ,[/ 7 ) of the Vlasov-Poisson 
system with the following properties: 

(i) / 7 (x,f)=0(£)^(7P), x,t>GlR 3 . 

(ii) (f°,P°,U ) = (f( h po,Uo), (/ 7 ,p 7 ,?7 7 ) is axially symmetric for |tI < To, 
more precisely, for all x,v G IR 3 and AgS we have 

P{Ax,Av) = P{x,v), p 7 (Ar)=p 7 (x), W(Ax) = U\x), 

and (/ 7 ,p 7 ,f/ 7 ) is not spherically symmetric for 7 7^ , i. e., the above 
identities fail if S is replaced by 5*0(3). 

(iii) p 7 G C^IR 3 ) and £/ 7 GC b 2 (IR 3 ) for | T | <To- 

(iv) The mappings}— 7o,7o[3 7 ^ P 7 an d ] — 7o>7o[3 7 l— ¥ U 1 are continuous 
with respect to the norms \\-\\ loo or 1 1 • 1 1 2 00 r &spectively. 
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The question whether the axially symmetric steady states obtained above 
are static or not is addressed at the end of this section. In order to prove the 
above theorem we first deduce the semilinear elliptic problem ( |1.2p introduced 
in the introduction: 

Lemma 2.1 Let f be given by f(x,v) :=4>(E)if)( , yP) for some potential U : 
1R 3 -^]R. Then 

p(x) = h(j,r(x),U(x)), xGlR 3 , 

where 



r(x):=\Jxl + X2, xGE 

and 



h(j,r,u) 



rE rj2(E-u) 

2tt / <p(E) / ^rs)dsdE , u<E 0l 7GH, r>0, 

Ju J-y/2(E-u) 



, m>4 7eM, r>0. 



Moreover, h G C^R x [0,oo[xIR), d r /i G C 1 (IR x [0,oo[xK) ; and for every 
bounded subset B C IR x [0,oo[xIR i/iere exzsi constants C>0 and ue]Q,1[ 
with 

|9 r /i(7,r,w)| < Cr, 
|/i(7,r,«) — /i(7',r,M')| < C(|7 — 7'|r + |u — 
|3 u A( 7 ,r,iO-d u /i(7Vy)| <C(|7-7'| + |«- W r) 
for all ( r y,r,u), (Y,r,u')&B. 

Proof: The formula for h follows by introducing cylindrical coordinates with 
respect to (—xz,Xi,0)/r(x) in velocity space; if r(x) = then ipi^P) = 1 and 
one can use spherical coordinates. The function h is easily seen to be con- 
tinuously different iable with 



pE ry/2(E-u) 

djh("y,r,u) = 2nr I (f>(E) / sijj'('yrs)dsdE, 

Ju J-y/2(E-u) 
rE ry/2(E-u) 

d r h(^,r,u) = 27T7 / (f>(E) / sip'('yrs)dsdE, 



-y/2(E-u) 

Eo ( fWT^ u fWT^ 7\\ 0(E)dE 



d u h(y,r,u) = -2tt / {^{ ir p{E -u)) +^{- 1 r^2{E -u))) 



2(E-u) 



for u<Eq. The assumptions on ip imply that \ip'(P)\<C\P\ on bounded 
sets containing 0, which yields the estimate for d r h. The second estimate is 
straightforward. Since d u h is continuously differentiable with respect to 7, it 
is locally Lipschitz with respect to 7. As to the asserted Holder continuity of 
d u h with respect to u, take ( r y,r,u), ( , ~f,r,u')^B, B C IR x [0, 00 [xIR bounded, 
assume u<u' and let - + - = 1. Then 



\d u h(~f,r,u)-d u h(7,r,u)\<C / 

Ju 



u 

<*<> - . <f){E) 



+C { ° )y/E^u-y/E^v/\-^=(lE 

Ju' 

+C 



yfE-u 
4>{E)dE 



\JE — u y/E — u 

< C\u-u'\*~2 +C\u-u'\i + C\u-u'\v~^ 



I 

,, 1 ,,11 



where we have used Holder's inequality for the first and the last term. Since 
i — I > this completes the proof. • 

We note that the above estimates would simplify if we asssumed that <fi is 
Holder continuous, but this would exclude the polytropes with — |</x<0. 
Next we collect some further properties of the spherically symmetric steady 
state (f ,po,U ): 

Lemma 2.2 The spherically symmetric steady state (fo,Po,Uo) has the fol- 
lowing additional properties: 

(a) The function h(0,-,-) does not depend on the variable r(x); we will write 
it as h = h (u) for simplicity. 



(b) The potential U is given by 

U (x) = - I P^-rdy = ~ t A s 2 p Q (s)ds-4n f°° sp (s)ds, 16E 3 , 
J \x — y\ \x\Jo J\x\ 



and 

4% r\ x \ 



47r r\ x \ 

U' Q (\x\) = TT : / S 2 Po (s)ds, XGlR 3 . 



(c) p is decreasing with po(0)>0, Uq(0)>0, for every R>0 there exists 
C>0 such thatU' Q {r)>Cr, re[0,R], andU (l) = E . 
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(d) p' is Holder continuous, and UqEC 2 (M 3 ) where IR 3 := IR 3 \ {0}. 

We identify p and Uq as functions of \x\ with p and Uq as functions of x; 
the derivative with respect to \x\ is denoted by ' . 

Proof: The assertion in (a) is obvious from Lemma |2.1| . Since we require 
that Ywi\ x \^ 00 Uq{x) = 0, the assertion in (b) holds by uniqueness. Since ho is 
decreasing and Uo is increasing we find that po is a ls° decreasing, and since the 
steady state (fo,Po,U ) is assumed to be nontrivial we must have po(0) >0. 
Thus actually U' Q {r) > 0, r > 0, and since Uq(0) = ^Po(O) > this implies the 
estimate on Uq from below. The assertion that Uo{1) = Eq follows from the 
form of ho and the fact that by assumption supppo = -Bi. The regularity of 
Uq follows from the formula above and the fact that poGC^(IR 3 ). Finally, 
the Holder continuity of p' Q (r) = d u h (Uo(r)) Uq{t) follows from Lemma P7T| . 

• 

We want to find solutions of the equation 

AU = 4irh(7,r(x),U), (2.1) 

and the central idea is to reformulate this as a problem of finding zeros of 
an operator T which acts not on the space of potentials directly but on 
deformations of the given spherically symmetric potential Uq. We now define 
the Banach spaces which will serve as domain and range of T: 

X:={feC s (B 3 )\f(0) = 0, feC\B 3 ), 3C>0:|V/(x)|<C, xeB 3 , 

Vx G Si : lim V f{tx) =: V/(0rr) exists, uniformly in x G Si \, 



which we equip with the norm 



x 



■ sup \Vf{x) 



/ex, 



and 



F:={/eC s (B 3 )|/(0)=0, f£C l (B,), BC>0: |V/(x)| <C|i|, x£B 3 , 

Vl€Sl : lim YlM = ..YI<M exists , uniformlyinieSl i 

t\o t J 
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which we equip with the norm 



x£B 3 



\X\ 



For / e X the function V/(0-), being defined as the uniform limit of functions 
in C(Si), is itself in C(S 1 ). Furthermore, since /(0) = 0, 

\f(x)\ <\\f\\ x \x\, xeB 3 , feX, 

and the norm || • \\ x is equivalent to the norm 



x-= sup 

X&B3 



\x\ 



■ + |V/(x)| + sup|V/(0x)|. 



It easily follows that (X, \ \ • \\ x ) is a Banach space. For / e Y note first that 

/M :=lim /M = IWM. X 

2 t\o t 2 2 



uniformly in xESi. We have 

v/(o-) /(o 







o 2 



eC(^), |/(x)|< 



|x| 2 , x<EB 3 , 



and the norm 



is equivalent to the norm 



y := sup 



\M\ , |v/(x) 



■ + 



+ sup 

a;GSi 



/(Ox) 



2 



+ 



V/(0x) 







from which it follows that (Y, \\ ■ \\ Y ) is a Banach space. 

Using the elements in the Banach space X we can deform spherically 
symmetric sets, e. g., the level sets of the given, spherically symmetric static 
solution, into axially symmetric sets in the following way: 



Lemma 2.3 For(eX define 



x 



g^.B^lR 6 , g c (x): = x + C(x) — , x€B 3 , g c (0):=0. 



Then there exists r > such that for all 

(eQ-.= {(eX\ 

the following holds: 



<r} 



10 



(a) 0£ : B 3 — > B 3 ^ : = g^(B 3 ) is a homeomorphism, g^: B 3 ^ B 3 ^ is a C x -diff- 
eomorphism whose Jacobian satisfies the estimate 

\Dg ( (x)-id\<^, xEB 3 , 

and for every x E Si the restriction 

g c : 0,3x 3y^g c (y)E 0, \g((3x) k 

zs one-to-one, onto, and preserves the natural ordering of points on the 
line segment 0, 3x. 

1 3 

(b) ^\x\<\g c (x)\<-\x\, xEB 3 , and g c (B 1 )cB 2 , B 2 Cg c (B 3 ) C B 4 . 

(c) g c (Ax)=Ag c (x), xeB 3 , and g^ 1 (Ax) = Ag^ 1 (x) , xEB u , AeS. 

(d) \Dg^ 1 (x) — id\ <-, xEB 3 ^, and there exists a constant C>0 such that 
for allCX'ett, 

^\g c (x)-g c ,(x)\ + \Dg c (x)-Dg c (x)\<C\\(-C\\ x , x E B 3 , 

and 

\g^(x)-g^(x)\<C\\C-C'\\ x \x\, xEB 2 . 
Proof: On B 3 we have for i,j = 1,2,3, 

d Xi g u (x) = 5 tJ + d Xi ((x) ^ + <M ^ _ ^ (2.2) 

which implies that 

\Dg c (x)-id\<3\\(\\ x <^, xEB 3 , 

provided < 1/6. Using the inverse function theorem we obtain the first 
two assertions in (a). Since g^(y) G]0,oo[y for every yEB 3 the remaining 
assertion in (a) follows. The assertions in (b) are obvious, provided r > is 
chosen sufficiently small, and so are the assertions in (c). As to (d), the first 
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assertion follows by choosing r still smaller, since Dg^(x) = (Dg^)" 1 ^^ (x)) 
The estimate for g^ — g^ is immediate from the definition of g$. The estimate 
for Dg^ — Dg^i follows from Q2.2| ). Finally, ieB 2 implies that x£g^(B 3 )r\ 
g^(B 3 ), and there exists yEB 3 such that x = g^(y). Thus 

I^V)-^ 1 ^)! = \9{ 1 (gc(v))-y\ = \9{ 1 (g<;'(y))-g{ l (gdy))\ 

<2\g c (y)-g c (y)\<2\\(-C\\ x \y\<4\\(-C\\ x \x\ 

by the mean value theorem, the estimate for Dg£ which we already estab- 
lished and the fact that g^(y), gc,'{y) Cg^(B 3 ). • 
We want to find solutions of the reduced problem (|2.1|) of the form 

U(x) = U c (x):=U (g^(x)), x6B 3 ,(, 

for some C^^- Of course U will have to be defined on all of IR 3 , but this 
will be easy once we have it on -83,^. Using the fundamental solution of the 
Poisson equation we integrate fl2.1|) and transform our problem to that of 
solving the equation 

TTfs.f h ( r f^(y)Mg^ l (y))) 

U {x)+ : — — — dy = 0, x<EB 3 ; 

jb 3,< \g(\ x ) y\ 



observe that is invertible. It turns out that we can avoid the dependence 
of the domain of integration on (, and also that the operator above is not 
quite the right thing yet. We are now in the position to give the proof of the 
theorem: 

Proof of the Theorem: For ( G f2 and 7 G IR we define 

rrt r\( \ tt ( \^( h (l^(y),U (g^ 1 (y))) 

T (7,C) 0) := U (x) + / — ^ — ^ dy 

Jb 2 \g c (x)-y\ 

-U o {0)- — ^ dy,x<EB 3 . (2.3) 

Jb 2 \y\ 

Assume we already know that this defines a continuous operator T :] — 
l,l[xO— ►V, that T is continuously Frechet differentiable with respect to 
C, and that 

<9 ( T(0,0):X^r 
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is an isomorphism; that all this is indeed the case is shown in Sections 3 



and 4, cf. Proposition |3J] and Proposition |4.1| . Note that by definition of Y 
we must have T(7,£)(0) = which is the reason for subtracting the constant 
term in the definition of T. By assumption (02) we know that T(0,0) = 0; 
note that go = id and that suppp = supph Q oU = BxC B 2 . By the implicit 
function theorem there exists a constant 7qG]0, 1[ and a continuous mapping 

such that 

T(7,C 7 )=0, 7e]-7o,7o[ 
and C° = 0, cf. [§, Theorem 15.1]. Let C = C 7 f° r some je] — 7o,7o[ an d define 

p c (x):^h^,r(x),U (g^(x))), xEB 2 . (2.4) 

Then p^ECs{B 2 )\ at the moment the differentiability at x = is not yet 
obvious. By Lemma ^]T]p^(a;)>0 if and only if Uo(g'^ 1 (x)) <Eq which by 
Lemma |2.2| is equivalent to |(7^ 1 (x)| < 1. Therefore, by Lemma [273 

o 

snppp ( = g ( (B 1 )cB 2 , 
and we can extend p^ by to all of IR 3 , obtaining an element of C C (IR 3 ) with 

o 

suppp^Ci?2 which for the moment need not satisfy (|2.4|) everywhere. The 
equation T(j,Q = can now be written as 

Uo(x) = -[ dy+C,xEB 3 , 
Jb 2 \g c {x)-y\ 



or 



where 



Now define 



U (g^(x)) = -f £MLdy+C,xeB 3> 
s jb 2 \x — y\ 

C:=l/ (0)+/ P -Mdy. 
Jbo. \y\ 



U c (x):=- f P^jdy + C, xEJR 3 
^ v ' Jn 3 \x-y\ 
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Clearly, L^gC^H 3 ) with 

U ( (x) = U (g^(x)), xeB 2 cB u . (2.5) 

This immediately implies that p c G C^IR 3 ) and U{ G C& (R 3 ) with A?7 C = 47rp c 
on IR 3 . On the other hand, 

A£7 f (z)=47r%,r(aO,E/ c (aO), ^B 2 C% (2.6) 

If fact the latter equation holds on all of IR 3 . To see this we have to show 
that p^(x)—h(j,r(x),U^(x)), zGlR 3 , that is, U^{x) > E for x G IR 3 \g f (Bi). 
We know that 

AU c (x) = 0, xGlR 3 \g f (B 1 ), 

and 

lim U c (x)=C, U c {x)=E , xedg c (Bt), U c (x)> E , xeB 2 \g c (B 1 ), 

\x\— >oo 

where we used the identity (|2.5|) and the fact that U is strictly increasing as 
a function of \x\ with Uq(1) — E q . The assumption C <E would contradict 
the maximum principle. Thus, C > E , and again by the maximum principle, 
U ( > E Q on IR 3 \g c (Bi). Therefore, fl2J) does hold on all of IR 3 . 
If we define p 7 '■ = p^, U 1 := Uq , and 

P(x,v) : = 4>(^v 2 + U' Y (x)^ij(j(x 1 V2-X2V 1 )), x,velR 3 , 

then the assertions (i)-(iii) of the theorem are established, except for the 
assertion that the solution is not spherically symmetric for 77^0. To see 
the latter, choose a;GlR 3 with p 7 (x)>0, x%^0, x 2 = X3 = 0. There must 
then exist r/^0 such that \rf + U~ 1 (x) <E . Let v— (0,0, r/) and v' — (0,77,0). 
Then there exists a rotation A around the xi-axis such that Av — v', and 
clearly, Ax = x. But since E(x,v ) = E(x,v ') and P(x,v) = 0^Xirj = P(x,v') 
we have 

P(x,v) =<p(E(x,v))i)(~iP(x,v))=<p(E(x,v)) = <p(E(x,v')) 
^4>(E(x,v'))^( 1 P(x,v')) = P(x,v'), 
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provided 7 7^0. The continuity properties asserted in (iv) follow from the 
fact that C 7 depends continuously on 7 with respect to the norm || • First 
the estimate 

\U^x)-U^\x)\<\\Ul ) \\j9^\x)-g^)\<C\\C 1 -C Y \\ x , xeB 2 , 

and the relation p~'(x) = h{j ,r(x) , £/ 7 (x)) imply that p 7 depends continuously 
on 7 with respect to || • H^. Since 

lP(x) = -f P^dy + U (0)+[ ?-iy±dy,xeJR 3 , 
Jb 2 \x — y\ jb 2 \y\ 

this implies that U 1 depends continuously on 7 with respect to || • \\ x OC) . Dif- 
ferentiating the above formula for p 7 we obtain the asserted continuity of p 7 
with respect to || • 1^ ^ and thus also of U"' with respect to || • || 2oo , and the 
proof of the theorem is complete. • 
Remark: 

(a) For fixed ip the family (/ 7 ,p 7 , U 1 ) is in a neighbourhood of 7 = unique. 
However, different functions ip give different families of stationary so- 
lutions in which (fo,po,Uo) is embedded. 

(b) The mass current density is given by 

p{x) := J vf y (x,v)dv 

rE rj2{E-W{x)) 

= 2tt / <f>(E) / , siphr(x)s)dsdEe t (x), x G IR 3 , 

JU->(x) J-y/2(E-W{x)) 

where 

e t(x) :=— -(-x 2 ,xi,0) 
r(x) 

denotes the unit vector field tangent to the orbits of points under coun- 
terclockwise rotations around the a^-axis; note that the integral above 
vanishes on the X3-axis, where e t is not defined. Now the average ve- 
locity j 7 /p 7 of the steady state vanishes identically if ip is an even 
function, in which case we obtain a family of static, axially symmet- 
ric solutions. In general, for example if ip(—P) < 1 <ip(P), P> 0, the 
average velocity does not vanish, and we obtain a stationary stellar 
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system which rotates around the X3-axis. It is also easy to see that 
the average velocity vanishes at the boundary of the support of p 7 . As 
already mentioned in the introduction, static solutions of the Euler- 
Poisson system are necessarily spherically symmetric. As we now see, 
a corresponding result does not hold for the Vlasov-Poisson system. 

(c) It would require very little additional efford to show that T is continu- 
ously different iable also with respect to 7. This would imply that the 
family (/ 7 ,p 7 ,£7 7 ) depends on 7 in a differentiable way with respect to 
the appropriate norms, cf. Cor. 15.1]. 

(d) If in our ansatz / depends only on the particle energy E then the right 
hand side of ( |1.2| ) does not explicitly depend on x. One can then apply 
a result by Gidas, Nl, and Nirenberg, cf. f| Theorem 4], to conclude 
that under mild regularity assumptions a corresponding steady state 
with finite radius and finite mass must be spherically symmetric with 
respect to some point in IR 3 . Therefore, it is necessary to include 
further invariants in the ansatz in order to obtain stationary models 
which are not spherically symmetric. 



3 The Frechet-difFerentiability of T 

The aim of this section is to prove the following result: 

Proposition 3.1 The mapping T :} — l,l[xQ — >F defined in ^2.3j is contin- 
uous and continuously Frechet- differentiable with respect to ( with Frechet 
derivative 



[«9<T(7,C)£](*) 



d u h( 1 Xy)My))vudy)--^^agt 1 (y))dy 



Z2\\g((x)-y\ \y\J u ' ' I^G/)! 



'B2\g((x)-y\ 3 ' ' |x| 
where 7 e]-l,l[, C^fi, (el, andU ( (y):=U (g^\y)), yeB 2 . 
In order to prove this result we need more information on the elements of 



the space X and the deformation mappings constructed in Lemma 2.3 
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Lemma 3.2 Let C^O. Then the following holds: 

(a) \((x)-ax')\<\\(\\ x \x-x'\, x,x'eB 3 . 

(b) ForxESi the mapping [0,3] 3t\— >C(t x ) ^ s continuously differentiable, 
and lim t \o — VC(0x) • x, uniformly in xESi. 

(c) The mapping [0,3] x S\ 3 {t,x) i— > V((tx) is uniformly continuous. 

(d) For xGSi the limits lim f \ =: and lim t \o D g^(tx) =: 
Dg^(Ox) exist, uniformly inx(zS\. 

Proof: The assertion in (a) follows easily by distinguishing the cases G x, x' 
and O^z, i'. As to (b), 

—£{tx) = VC(tx) • x -> VC(0a;) • x, t \ 0, 

by definition of the space X, and the rest follows. The assertion in (c) follows 
from the fact that VC eC(5 3 ) and VC(te) — ► VC(0x) uniformly in x G Si as 
t\0. The assertions in (d) are easy consequences of the definitions of 
and the space X, together with (|2.2|). • 

Next we establish some estimates for the spatial density induced by a defor- 
mation of the potential Uq: 

Lemma 3.3 For 7 g] — 1,1[ and ( GO let 

Py,d x ) '■= h {.li r i x )iUa(g^ l (x))), xeB 2 . 

Then the following holds: 

o 

(a) p 7j £ G Cs(B 2 ) PlC (B 2 ) with suppp 7i ^ C B 2 , and there exists a constant 
C > such that for all 7 G] — 1 , 1 [ and ( G O, 

|Vp 7) ^(x)| <C|x|, xEB 2 . 

(b) There exists a constant C>0 such that for all 7,7' G] — 1, 1[ and(,('^ 

\py,d x )-py,c'( x )\<c(\i-i'\ + \\C-C\\ x )\ x \^ ^5 2 . 
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Proof: Lemma and Lemma [2.1| imply that p = p 1 ^ £ ^5(^2) nC^-E^)- 
For x G B2 we have 

Vp 7iC (x) = d r h(-f,r(x),U (g^ 1 (x)))'Vr(x) 

+ d u h( 1 ,r(x),U (g^(x)))VU (g^\x))-Dg^\x), 

and Lemma |2j|, the fact that U G C 2 (IR 3 ) with VU (0) = 0, and Lemma ^3 



imply the estimate 

|Vp(x)| <C\x\ +C\VUo(g^(x))\ <C\x\ +C\g^(x)\ <C\x\, xeB 2 - 

note that the range of U is bounded. Since x ^g^(Bi) implies U (g^ (x)) > 
E and thus p(x) = 0, the assertion on the support of p follows by 
Lemma ^]3| (b). The assertion in (b) is immediate from Lemma ^73] and 
Lemma |2.3| (d). • 
We shall need the following assertions on Newtonian potentials: 

Lemma 3.4 Let o-eCs(B 2 ) be such that 



\<J(X) 



■ sup ■ 

x£B 2 



< OO 



\X\ 



and define 



VJx) 



dy, xeJR 3 



Ib 2 \x-y\ 

Then V a G C 1 (1R 3 ), and there exists C>0 such that for all a as above the 
following estimates hold: 

(a) \W a (x)\ <Cc a \x\, xGlR 3 , 

(b) \VV a (g c (x))-VV„(g^x))\<Cc„\\(-C\\¥\x\, xeB 3 , CC'eH. 
Proof: For aeC s {B 2 ) we have VK(0) = and thus 



|W,(z)|< 



x-y y 



iB 2 \\x-y\ 3 1 2/ 1 3 
Let Xy^O and r:=2|x|. We obtain the estimate 



a(y)dy 



x g nr. 



|VK(x)| < c a J^ 



B 2 \B r 
■■h+h- 



x-y 



\x-y\ 3 \y\ 3 



\y\dy + c a 



' +T7At)\v\dV 



B 2 nB r \\x-y\ 2 \y\ 
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For almost every y G B 2 there exists r G [0, 1] such that 

4 



x-y + y 



\x-y\ 3 1 2/ 1 3 



< \x\ 



rx — y\- 



and since for \y \ >r, 



r \y\ 

\TX-y\>\y\-\x\ = \y\-->Y> 



we can estimate the first term as 

h<Cc a \x\( -\-dy = Cc a \x\; 
Jb 2 \y\ 2 

constants denoted by C may change their value from line to line or even 
within one and the same line. For the second term we have 

h<1cA [ | — ^—r^dy+f r^cfo/) <4c CT / -\rdy = Cc CT r = Cc CT \x\, 
\JB r \x-y\ z JB r \y\ z J JB r \y\ z 



and the proof of part (a) is complete. As to (b) we have 



\VV a (g c (x))-VV a (g c ,(x))\<c a 



B 2 



\gd x )-y\ 3 \gc( x )-y\ 3 



\y\dy. 



Let xEB 3 and 5:= ||C — C'llx < 1> r i • = 2<5 , and r 2 :=4|x| >r\; recall that 
we chose the radius of the set f2 less than 1/3. We split the integral above 
into three parts, I±, I 2 , and I 3 , according to the decomposition 



B 2 = (B 2 \B r2 )u ((B 2 HB r2 )\ B n (g ( (x))) U (b 2 n B n (g c (o 
As to Ii we find for almost every yeB 2 ar between ((x) and ('(x) suc h that 



gd x )-y 9c( x )-y 



< 



c 



\ x + T Tkr y \ 



\g c (x)-y\ 3 \g^{x)-y\ 3 
note that both g^(x) and g?(x) lie on the line IRx. Since 

\((x)-('(x)\<\\(-('\\ x \x\=6\x\, 



\C(x)-C(x)\; 
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and for \y \ >r 2 , 



X 




x + t— — y 




\x\ 





X 



y-g c (x) + (C(x)-r) — 



>\y\-\gdx)\-\C(x)-C(x)\ 



3 5 5 3 

> \y\ - 2 M ~ 5 \ x \ ^ \y\ - 2 N = M ~ g r2 ^ § M> 



we find the estimate 



/i<c<*,||C-C'y*l 1^=^11^-^11x1*1 

To estimate the second term I 2 we start like for Ii, but for y^B ri (g^(x)) 
obtain the estimate 



X 






x + t— — y 




y 


\x\ 







X 

I x I 



>\y-gdx)\-\C(x)-C(x)\ 



> \y-gd x )\- 6 \ x \ ^ ij\y-9c( x )\- 



On the other hand for y 6 £? r2 we have 



3 3 



and 



J 2 < CcJ\x\J^ 



1 2r 2 
- — — -dy = CcJ\x\Anhi — 

B 2r2 (g ( (x))\B ri (g c (x)) \g C {x)-y\ 6 T X 



Cc a 5\x\\n-<CcJ 1/2 \x\ = Cc (T \\(-C\\ 
o 



/ 1 1 1/2 | 
X 



x\. 



As to the third term we have 



h < 2c n 



dy 



b. 



<4c ff / 
Je 



lOjcW) \gd x )-y\ 2 

dy 



+ 



By 



,(9<W) I^C'(^)-i/| 2 / 
Cc^ri = Cc (7 5|:r| =Cc (J ||C-C / || x |x| 



and the proof of part (b) is complete. • 
We are now ready to prove part of the assertion in Proposition |3.1| , namely: 
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Assertion 1: For 7 g] — 1,1[ and C G Q we have T(7,C) £ an( ^ ^ e mapping 
T:] — 1,1 [xQ — >Y is continuous. 
Proof: Let 

Vrt(x):=-f ^\dy,xem 3 , ( 7 ,c)e]-i,i[xn. 

Jb 2 \x — y\ 

The assertions in Lemma pO| (a) imply that G C 2 (IR 3 ) with VV 7i< j(0) = 0. 
Since 

T( 1 ,C)(x)=U (x)-V yX (g <: (x))-U (0) + V^(0), xEB 3 , 

cf. O, we have T( 7 ,C)(0) = and T( 7 ,C) G C 1 ^) nC 5 (fi 3 )- While we 
show that T(7,£)gF for (7,^) G] — l,l[xfi the arguments 7 and £ remain 
fixed, and we write V = Vy X . From 

VT("f,()(x)=VU (x)-VV(g ( (x))Dg c (x), xeB 3 , 

we obtain the estimate 

\VT( 1 X)(x)\<\\D 2 U \\Jx\ + 2\\D 2 V\\Jg ( (x)\<C\x\ 

with some constant C which depends on Uq and V but not on x. In particular, 
this shows that T( 7 ,£) G C 1 (i?3). Now fix x G Si. Since any point on the line 
segment 0,g^(tx) can be written in the form g^(rx) with rG [0,t] we have 

<9 x ,T( 7 ,C)(te) a^Z/oCte) 1. 



t t 
d Xi U Q {tx) 1 



DV(j ( (n)) J( (te))-9 IiJ( (te 



V^C/ o (0)-a;- I £>V(0)^M I -d Xi g^x) 



as t \ 0, uniformly in x G Si, by Lemma [T2] (d). This shows that T(7,C) G F. 

To show that T is continuous we fix ( 7 ',£') e] — 1, l[xfi. Constants de- 
noted by C may depend on (7', CO but not on (7,C) s] — 1, l[xfi or x&B 3 . 
Restoring the subscript of V we have 

||T( 7 ,C)-T(7 , ,C0lly = sup ^ VV lX (g c {x))Dg c (x) - VV jf>c ig c {x))Dg c (x) 



x£B 3 



\X\ 



< sup --(h + h + h), 



X 
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where for x G B%, 

h 
h 
h 



\Dg c (x)\\VV 7)C (g c (x))-VV y>c (g c (x))l 
\Dg c (x)\\VV y , c (g c (x))-VVy, c (g c (x): 
\VVw(gc(x))\\Dg € (x)-Dgc(x)\. 



Using Lemma |3.3| (b) and Lemma |3.4| (a) with cr = p 7 ^ — py £ we find 

|Vy 7 , c (^(x))-Vy v , ( ,(^W)l<C , (l7-7 , | + IK-C'llx)k(^)l 
and thus by Lemma |2.3|, 



Since Fy iC / G C 2 (IR 3 ) with Wy,C'(0) = we have by Lemma ^ (d), 
l2<C\g c (x)-g c ,(x)\<C\\(-C\\ x \x\, xEB 3 , 

and 

h<\\DX>^\Jg<>(x)\\\C-C\\x<c\\C-C\\ x \xl ^b 3 , 

and the proof of Assertion 1 is complete. • 

To deal with the differentiability of T we have to investigate the integrand 
of the first term in the formula for d^T, cf. Proposition |3.1| : 

Lemma 3.5 For 76] — 1,1 [ ; C^^; an d (Gl define 



o-^,d x ) ■=d u h(-i,r(x),U c (x))VU Ci {x) 



9c 



^g c \x)), xeB 2 , 



where we recall that U^(x) = Uo(g ( ~ (x)), x<EB 2 . Then a 1) ^^£Cs{B 2 ), and 
there exists C > such that for every 7 g] — 1, 1[ ; C G fi, and (GX, 

\<Ty£,s(x)\<C\\£\\ x \xl xeB 2 . 

Moreover, if we fix (j',C) G] — l,l[xfi there exists for each e>0 a 5>0 such 
that for all (7,0 G] -l,l[xfi with \l-l'\ + \\(-C\\ x <6 > and^eX, 



\°i,u{ x )- a i'£,tit x )\< e \\£,\\x\ x l X ^ B 2- 
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Proof: Since the range of Uq and thus also of Uq is bounded, the first factor 
in er-y,^ is bounded, uniformly in 7 and (, and the same is clearly true for 
the second and third factor. Together with 



\tty\x))\ < U\\ x \g^(x)\ <2\\Z\\ x \x\, xeB 2 , 



(3.1) 



the estimate for c 7 f g follows. The continuity of a 7 £ £ on _B 2 is clear, and 
at x = it follows from the estimate above. The symmetry follows from 
the corresponding properties of U , g^, and £. In the following C denotes a 
constant which may depend on Uq and (7', CO but not on 7,C;£ ; or x. Making 
excessive use of the triangle inequality we find that 



a 



7.C.C 



X) —Gryl££{X 



<c 



9 u / i ( 7 ,r(x),^(a;))-9 u / i (7 , ,r(a;),^(a;)) \^\x))\ 

■civc/c^-v^^lie^c -1 ^))! 



4c 



m\x) 



=:/i + / 2 + /3 + /4, xeB 2 . 

Now the estimate flOD together with the properties of the function h stated 
in Lemma EO imply that 



h<C{\ 1 -i\^U,(g- c \x))-U^,\xW)\\i\\M 
<c{\i-i'\^gf{x)-g^{x)\ v )\\^ x \x\ 

^C^-yi + IIC-C'll^lieLN, ze5 2 . 

The crucial estimate is the one for J 2 : it is at this point that we need the 
limit condition in the definition of the Banach space X and its consequences. 
First note that 

h < C\VU Q {gf{x))-VU Q {g^{x))\ \\t\\ x \x\ 
+ C\Dg- c \x)-Dg- cl \x)\\\i\\ x \x\ 

<CK-C\\ x U\\ x \x\ + C\Dg- c \x)-Dg- c ,\x)\U\\ x \x\- 
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Now with z:—g^(x) and z' : = g,, (x), 

\Dg^(x)-Dg^(x)\ = \{Dg c )'\z) - (Dg c )-\z')\ 
<C\Dg c {z)-Dg c ,{z')\ 

<C\Dg c (z)-Dg c (z)\+C\Dg C r(z)-Dg c (z') 
<C\\(-('\\ x + C\Dg c ,(z)-Dg c ,{z')\, 



and it remains to estimate the last term in the line above. From ( |2.2j ) we get 
the estimate 



\Dg c (z)-Dg c (z')\ < C\Vf(z) - VCV 

+ C|VCV 

+ C|CV)I 
=: Ji + J 2 + J 3 + Ja- 



z 


z' 


\z\ 


~\A 


1 


1 


\z\ 





+ £|C'(*)-CV) 

At 



+ max 

ij'=l,2,3 



ZiZj Z i Z j 



|z| 3 Iz'l 3 



With x : = x/|x| there exist s,s' > such that z = g^ (x) = sx and z' = g c } 
s'x so that s = kl, s' = k'l, and 



\s — s 



1 = ||z|- < -g^(x)\ <C\\C-C'\\ X , xeB 2 . 



Now given e > we can choose 5 > according to Lemma |37^ (c) such that 

IIC-C , llx< 5im P lies 

J 1 = C\VC(sx)-VC(s'x)\<e, xeB 2 . 
Using Lemma |2.3| (d) and Lemma |3.2| we obtain 



J 2 <c^+^|^-^<^||C-C , llxN<c'llC-C / ll 



X ! 



J3<^\Z-Z>\<C\\(-C\\ 

\z\ 



X ' 



and 



^4<C|/| ( j^ + A, ) \z-z'\<C\\(-('\\ x , 



\z\ 2 \z'\ 2 
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so that finally 

l2<C\\(-('\\ x U\\ x \x\+CeU\\ x \x\, xEB 2 , 

provided ||C — C'llx^^- The remaining terms I 3 and J4 are much easier to 
estimate: 

<c\\c-C'\\ x U\\x\^ x ^b 2 , 

and 

and the proof of Lemma |3.5| is complete. • 
To continue with the proof of Proposition |3.1| we denote for fixed (7,()e 
] — l,l[xfi by L£ the right hand side of the definition of d^T(j,QCj £ &X. 
We now show: 

Assertion 2: L^C(X,Y) is a bounded, linear operator, and for all^EX, 

lim r( 7 ,c+tfl-r (7 ,0 

t^O t 

with respect to || • 

Proof: It is convenient to introduce the auxiliary Banach space 

Y:={feC s (B 3 )\f(0) = 0, feC\B 3 ), 

3C>0:|V/(x)|<C|x|, xeB 3 }, 

which we equip with the norm || • || y ; clearly, Y is a closed subspace of Y. 
Since we already know that T maps X into Y it is then sufficient to show 
that L G C(X,Y) and that the asserted convergence holds. To see the former 
define 

V ( (x):=- / -a^{y)dy, xGlR 3 , 

jb 2 \x — y\ 



and 

W(x):=- 

Ib 2 \x — y 



W(x):=-f -^—p 1)( (y)dy,xeJR 3 . 
Jb 9 \x — y\ 
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Then V^eC^Ht 3 ), WeC 2 (JR 3 ), and we can write 



x 



(LO(x) = Vt(g c (x)) - Vt(0) - VW(g c (x)) ■ —£(x), x e B 3 . 



x 



This implies that for (el we have L£eC\B 3 ), (L£)(0) = 0, and 



x 



(VL£)(x) = V^( 5c (x))^ c (x)-L>^( 5c (x))^ c (x)^e( 



x 



x\ 

-WW{g c {x))D ( ^ ) f (x) - Vlf^W) • V£(x), x G 4- 
Using Lemma 3J5 and Lemma |3]4] (a) we obtain the estimate 

\^L0(x)\<Cm x \9d^\+mx)\ + C\\D 2 W\\Jg c (x)\ (^ + 1 V£(s)l) 

<C||£U*|, xeB 3 . 
In particular, this implies that L£ is differentiable also at x = 0, and 

P*|| y <C||£|| x , £eX. 

The symmetry of L£ follows easily from the corresponding properties of V%, 
W, (, and £. In order to show that L£ is indeed the Gateaux derivative 
of T at (7,C) m the direction of £ we choose to>0 such that ( + t£,<EQ for 
|t| <to- Although this is in conflict with earlier notation it is convenient to 
abbreviate 

x 

9t(x)=gc+tz(x)=x + (((x)+tt i (x))—, xeB 3 , fe]-t ,*o[- 

Then for each 2G-B2 the mapping ] — to,to[3t\-^- g^ 1 (x) is continuously dif- 
ferentiable with 

^\x) = -{Dg t )-\g;\xmg-\x))f^\. 
dt \g t (x)| 

To see this, define for fixed zEB 2 the mapping 

G(t,x):=g t (x)-z, te]-t ,t [, xeB 3 . 
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Since G(t,g^ 1 (z)) = 0, t e] — £o,^o[, the asserted regularity of gf 1 with respect 
to t follows from the regularity of G, the fact that d x G(t,x) = Dg t (x) is invert- 
ible, and the implicit function theorem. If we now differentiate the identity 
x = g t (gi 1 (x)) with respect to t we obtain the formula for j|(? t -1 (:r). 
It will also be convenient to abbreviate 

p t (x): = p 1 ^ +t t.(x), a t (x):=G 1 £ +tU (x), te]-t ,t [, xeB 2 , 

and define 



F(t,x):= I ( Y ^—-^-\p t {y)dy, xeM 3 , te]-t ,t [. 
jb 2 \\x y\ \y\ I 



Then except for d\F all derivatives of F up to second order exist and are 
continuous on ] — t ,^o[xIR 3 , and 

d t F(t,x) = -f (-J—-l-\a t (y)dy, 
JB 2 \\x-y\ \y\J 

VF(t,x) = -f ^\ Pt {y)dy. 
Jb 2 I a: — 2/| 

o 

These results follow easily from the fact that ptECl{B 2 ) and 
f t Pt(y)=dMlAy)Mg^\y)))VU (g^\ y ))-j t g^(x) 

= -d u h( 1 ,r( y ),Uo(gr 1 (y)))VUo(gr\y))(Dg t )- 1 (g^(y)) 



! . A 9t '( .'/ ) 



tor(v))r^i 



\9t 



= -9 u M 7) r(y), f / (^ 1 (y)))V(^o(^ 1 (y)))-r L ^^^ 1 ^)) 

\9t (y)\ 

= -My)- 

Now 

r(7,C + ^)W-T( 7 ,C)W _ F(t,g t (x))-F(0,g,(x)) 

+ F(QM*))-F( M*)) ]te] _ toMxeB3] 
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and 



x 



m)(x)=d t F(0,g (x)) + VF(0,g o (x))-TlZ( x )> X ^ B ^ 



x 



one should be careful to note that here go = 5'c+o£ = <?<:• We claim that as 
t->0, 

F(t,g t (x))-F(Q,g t (x)) . xx . , 

v >wv n ^d t F{0,g (x)) (3.2) 



and 



F(0,g t (x))-F(0Mx)) ^ VF(0jflb(x)) . * m 
t \x\ 



(3.3) 

where both limits are understood with respect to the norm || • || y . This would 
then prove that L is the Gateaux differential of T at (7,0- As to ( |3.2|) we 
observe that 

v F(t '^ (a;)) ; F(0 '^ (a;)) -v(^(o,, (,)) 



VF(t,g t (x))-VF(0,g t (x)) 



Dg t (x)-\7d t F(0,g (x))Dg (x) 



< 



VF(^ t (i))-VF(0, ft (i)) 



-V0 t F(O,<fc(a;)) 



|£<7t(*) 



+ \Vd t F(0,g t (x)) - Vd t F(0,g o (x))\\Dg t (x) \ 
+ \Vd t F(0,g o {x))\ \Dg t (x) - Dg (x) \ 
=:h+I 2 + h. 

Let e > 0. For every z G IR 3 there exists r between and t such that 
VF(t,z)-VF(0,z) 



Vd t F{0,z) 



Vd t F(r,z)-Vd t F(0,z 
1 



b 2 \z-y\ 



vAy)-vo(y))dy 



and using Lemma [3]4] (a), the latter integral can be estimated byCe||£|| x |z|, 
provided 

Wr(y)-^o(y)\<^U\\ x \y\^ y£B 2 . 



Lemma |3.5| therefore implies that for 5 > sufficiently small we have 

h<Ce\g t (x)\<Ce\x\, xeB 3 , 
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provided \t\ <8. Note that C depends on ( and £, but not on t or x. Again 
by Lemma 3J3 and Lemma |3.4| (b) we find the estimate 

h<C\\i\\ x \\C + ti-C\\f\x\<C\t\ 1/2 \xl XGB 3 , 
and by Lemma |3.5| , Lemma |3]4] (a), and Lemma [2.3| (d) we conclude that 

l 3 <CU\\ x \go(x)\\\(+tt-(\\ x <C\t\\x\, xeB 3 . 



This proves the convergence in 
observe that for every x G B 3 , 

F{0,g t {x))-F{0,g o {x)) d 



with respect to 



As to (O) we 



x 



t dt 
where r lies between and t g] — £o?M- Therefore 



-F(0,g t (x)) lt=T = VF(0,g T (x)) ■ j-£(x) 



,F(0,g t (x))-H0,9o(x)) 



t 



V(VF(0,g (x))~a< 
v \x\ 



v[(VF(0,^(i))-VF(O l5 o(x)))~e^ 

<\D 2 F(0,g T (x))-D 2 F(0,g (x))\C\x\ 
+ \D 2 F(0,g Q (x))\\Dg T (x)-Dg (x)\C\x\ 



+ \VF(0,g T (x))-VF(0,g o {*))\ 



x 



x 



<C\t\\x\+C\D 2 F(0,g T (x))-D 2 F(0,go(x))\\x\, xe fig- 
Since D 2 F(0, •) is uniformly continuous on B4, which contains g T {x) for x G B 3 
and re] — t ,t Q [, cf. Lemma ^j3| (b), and 

\g T {x)-g {x)\<U\\ x \r\<C\t\, xeB 3 , 



we obtain the convergence in (|3.3| ) with respect to the norm || • || y . This 
completes the proof of Assertion 2. • 

Since a continuous Gateaux derivative is a Frechet derivative the proof of 
Proposition 3J. will be complete, once we show: 

Assertion 3: The mapping ] — 1, l[xf2 3 ~ * 9(T(jX) G C(X,Y) is con- 
tinuous. 
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Proof: Fix (7',^) e] - l,l[xfi and let (7,0 e] -l,l[xfi and ^el with 
||£|L = 1. Since 



[5 c r(7,C)e](x)-[9 c r(7',C , )e](x; 

Y 1 1 \ 



B2 [\\gt(x)-y\ \y\ 



^,c,dv)- 



.\gc( x )-y\ \y\. 



Ay) 



dy 



B 2 



\gd x )-y\' 



( x 9c( x )-y ( \ 



x 



\x\ 



we find 



V ([0 C T( 7 , C)e] (*) - [^(7', Of] (^)) I < Ix + J 2 + / 3 + / 4 + / 5 + / 6 



where 



/1 



1 



fi 2 VbcC^-yl M 
/ 1 



B2\\9((x)-y\ \g c (x)-y\ 

gd x )-y 



^1 ) ( <T 7,C,f(2/)- <7 'y,f'.f(3/)) d 2/ 

° 7 ',c>dy) d y 



:(p~f,dy)-PY,c(y)) d y 



B2 \gd x )-y\ 3 
( gd x )~y 9c( x )-y 



B2\\gd x )-y\ 3 \gt>( x )-y\ 3 , 
gd x )-y 



Py,c(y)dy 



\Z( x ) 



:(pi,dy)-pj',<:'(y)) d y 



B2 \gd x )-y\ 3 ' 
( gd x )~y 9c( x )-y 



x 



X 



B2\\gd x )-y\ 3 \9<'(?)-v\ l 



py',c(y) d y 



X 



\ \x\ , 



Given e > we choose 5 > so that the second estimate in Lemma |37| holds 
for I7 — 7' I + || £ — C'Hx < Then Lemma (a) implies, with z = g^(x), the 



estimate 



h< 



Defining 



1 



B2 \z-y\ 



°y,c,dy) - (T -y',c',dy)) d y D 9d x ) 



<Ce\z\<Ce\x\, xEB s . 



V(x) 



b 2 \x-y\ 
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we have V G C^IR 3 ) with W(0) = 0, and 



VV(g c (x))Dg c (x)-VV(g c (x))Dg c (x] 
<C\vV(g c (x))-VV(g c (x))\ + \VV(g ( ,(x)\\Dg ( (x)-Dg c (x)\ 

<c||c-C't /2 N+ci^^Hlc-C'llx^cilC-C'llfN, xeb 3 , 



where we have used the first estimate in Lemma |3.5| , Lemma [3.4| , and 
Lemma |2.3| . In order to estimate the remaining terms we define 

1 



V^(x) 
Then V^eC 2 (M 3 ) and 



b 2 \x-y\ 



D 2 V 7)C (g c (x))-D 2 V y><l (g ( (x)) \D g< (x)\ \t(x)\ 



< C \X | 1 1 p 7jC - Py iC , 1 1 J /16 1 1 p 7jC - Py , C ' 1 1 ^ 4 1 1 Vp 7)C - Vpy >C 



3/4, 



1 3/16 



J 3- 



^^(iT-yi + IIC-C'llx) 1*1. 
where we have used Lemma |3.3| and Lemma 1] . Next we have 

h=\DX'A9d^))Dgc(x)-DX'A9c^))Dg C (x)\\ax)\ 

<C\x\\DX',C<9d^)- D2 yY,C'(9c(x)\ + C\x\\Dg c (x)-Dg cl (x)\ 

^elxl+ClxlHC-C'll^ X ^ B ^ 

provided ||C — C'llx * s sman enough, where we have used the fact that D 2 Vy^r 
is uniformly continuous on £> 4 E3 g^ (x),g^'(x) and \g^(x) — g^>(x)\ < \\( — CWx- 
By Lemma [T3| (b) and Lemma |3]4] (a) for cr — p 1 ^ — py ^ we obtain, with 

z = gd x )> 



b 2 \z-y\ 



<c(|7-yi + IIC-C'll)kl, xeB 3 . 



Finally, 

i 6 <c|w 7 ^(2 C (z))-w 7 ^^ 

<C||C-C'LM> ^5 3 . 
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We have shown that for every e>0 there exists 5>0, depending on (7',C')> 
such that for all (7,0 e] - 1, l[xfi with \l-l'\ + \\(-('\\ x <5 and all ^eX 
with \\£\\ x = 1 we have 

||a c T(7,C)e-5 ? T( 7 , ,C / )elly<e, 



and the proof of Proposition 3.1 is complete. 



4 <9(T(0,0) is an isomorphism 

The aim of this section is to prove the following result: 

Proposition 4.1 The mapping d^T(0,0) :X-^Y is a linear isomorphism. 

Let us abbreviate L £ := c^T(0,0)£ for £gX. In order to prove the result 
above we rewrite Lo£: Observe first that go = id, and therefore U^, defined in 
Proposition |3.1| , coincides with the potential Uq of the spherically symmetric 
steady state we started with, if C = 0. In particular, Po(\x\) = h (Uo(\x\j) = 
h(0,r(x),U o (\x\) for xGlR 3 , and 

p[,(|x|) = ^(0,r(x),f/o(kl))^o(kl) = ^MO,r(x),[/o(x))Vf/o^)-^, xeIR 3 . 
Therefore, 

{L£){x) = -f (-J—-*\M\ y \)t(y)dy- [ ^% P o{\y\)d V ~Mx) 
JB 2 \\x — y\ \y\J JB 2 \x — y\- i \x\ 

= -U r (\x\)Z(x) - I (-J—-l-\ p ( i (\y\)£(y)dy, xeB 3 , £gX 
JB 2 \\x-y\ \y\ I 



Now let 



™ (I) ~-c5wUi;=sr w\) p ' M)mdy ' x€Bl - ?eCs(B3) - 

Then we can write 

(L £)(x) = -U^\x\)[(id-K)£}(x), xeB 3 , £eX. (4.1) 



As a first step towards proving Proposition |4.1| we show: 
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Assertion 1: The linear operator K :Cs(B 3 ) ^Cs(B 3 ) is compact, where 
Cs(B 3 ) is equipped with the supremum norm \\ ■ ] \ oo . 

Proof: For £eC s (B 3 ) let 



VAx):=- 



b 2 \x-y\ 



:Po(\y\ny)dy,xeiR 3 



(4.2) 



Then V^eC l (M 3 ), W € (0) = 0, and 

1 



l%(\x\) 



mix) -mo)), xeb 3 . 



Using Lemma |2.2| (c) we obtain the estimate 

\(KO(x)\ < J-H|w € ||J*| <C||eiL, xeB 3 , 

where the constant C depends on p and U , but not on £ or x. Thus K maps 
bounded sets into bounded sets. We claim that K£ is Holder continuous with 
exponent 1/2, uniformly on bounded sets in Cs{B 3 ). Let M>0 and assume 
HClloo — M. In the following, constants denoted by C depend on p Q , Uq, and 
M, but not on £ itself. There exists a constant C > such that 



|V^(x)-V^(xO|<q|p^||Jx-xf/ 2 , x,x'eB 3 , 
cf. [g Probl. 4.8]. Since W e (0) = 0, Q implies 

\VVi.(x)\<C\x\ 1/2 , xeB 3 . 
Now let x,x' G-B3 and \x\ < \x'\. Then 



(4.3) 



\m)(x)-(KO(x')\< 



\Vt(x)-Vt{0)\ 



ir (\x\) u^x'i) 



and we obtain for some zEB 3 with \z\ < \x'\ the estimates 



h< C 



\U' {\x 


)- 


-U' {\x'\)\ 




x\ 


\x' 





X 



<C\x-x'\ l l\ 
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and 

so that 
Also, 



h<^\^V^z)\\x-x'\<^-\z\ l/2 \x-x'\<C\x~x'\ 1/2 



X 



\(K£){x) - (Kg){a/)\ <C\x-x'\ 1/2 , x,x' E B 3 . 



| {K£) (x) | < C\ VVs(z) | < C\x\ l ' 2 , x E B 3 , 

and we have shown that K maps bounded subsets of Cs{B 3 ) into bounded 
and equicontinuous subsets of Cs{B^). Thus K is compact by the Arzela- 
Ascoli theorem, and the proof of Assertion 1 is complete. • 

As second step in the proof of Proposition [4. 1| we show: 

Assertion 2: id — K :Cs{B 3 ) — > Cs{B 3 ) is one-to-one and onto. 

Proof: Since K is compact is suffices to show that id — K is one-to-one. Let 
£ eCs{B 3 ) with £ — K£ = 0. In order to show that ( = we expand £ into 
spherical harmonics Yi m , ZgINq, m = —/,...,/, where we use the notation of 
[p] , Ch. 3] concerning the latter. Denote by (r,9,(f)) and (s,t,i[)) the polar 
coordinates of a point x or yEB 3 respectively. For £gINq and m = —l,...,l 
we define 

T J\x\=r 

Using the expansion 

i^pE E ^lirjr.^M), 

l X i/l I=0m=-i ZI + ir > 

where r< : = min(r,s) and r> : = max(r,s), cf. |ll], Eqn. (3.70)], we find that 
4tt 1 



2Z + l^(r)Jo ruv 'VH> +1 s l+1 JJ\y\=s 
4tt 1 z-3 



2Z + l£#(r) Jo 



i 1 V 


-y| 








' > 










s l+l 



p'o( s ) I h+t ~~ ) / , Y L{ T ^)i{y) dS v ds 
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This implies that 

too(r) = — 777TT p'o(s)s(s-r)£ 00 (s)ds. 
rUQ\r) Jo 

Clearly, £ o vanishes in the limit at r = 0. Let R > be such that £00 vanishes 
on [0,R). Then for re [R, 3], 

47r r 

\^o(r)\<— t^HPoIL su P I&oOOI / s (r-s)ds<C(r-R) sup |£ o(s)|- 



This implies that £00 vanishes in a right neighborhood of R and thus on the 
whole interval [0,3]. Up to multiplicative constants the spherical harmonics 
for / = 1 are given by sm.9e ±l< ^ and cos#, and the fact that £ G Cs implies that 
fi_i = fio = fii = 0. LeU>2. Then 

4-7T 1 / f r /"3 j-' \ 

6m(r) = - 2 ^ + 1£/ , (/ p' (s)—s% m (s)ds + J r p' ( s )— s % m ( s )ds\, 



and 



J-i ,3 J.-1 



^ * 27TT ^ ll6 "»- M ( -' ;)(s) ^ s3<fa+r i (-rt)W^* 

^ll&™IL(^/Vrt)(s)s 3 * + r^Vpi)(s)^) 

/ 1 3 r ^ 

■||6m|L( — ^ r3 Po( r ) + — / Po(s)s 2 ds + rp (r)ds 



'0\ 

4vr 1 
" 2Wfr) 

4vr 1 
~ 2/ + l^(r) 

= ^Yll6mlL, 

and since 2/ + l>3 this implies that £/ m vanishes for l>2 as well. We have 
shown that id — K is one-to-one as claimed, and Assertion 2 is therefore 
established. • 
It is clear that L : A — ► y is now one-to-one as well: just observe ( |4.1| ) and 
the fact that U' Q {r) > for r > 0. It remains to show: 
Assertion 3: L : A — > V onto. 

Proof: Let gGF and define q:=g/UQ. We claim that gG A. To see this we 
first observe that q EC 1 (B 3 )nCs(B s ), and 



Co(M) 2 M 
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By definition of Y and since Uq G C 2 ([0,oo[) with Uq(0) > we have that for 
every xES\, 



Vq(tx) 



Vg{tx) t g(tx) 



t U' {t) t 2 
Vg(0x) 1 g(0x 



x 



ttf(O) 







USi0) usW x 



as t \ 0, uniformly in x G S\. 

Since X C Cs{B s ) there exists by Assertion 2 an element £ G Cs{B 3 ) such 
that 

This implies that L ^ = g and thus that L is onto, provided (el. To see 
the latter we observe that £ = + g is Holder continuous since is Holder 
continuous. As above we conclude that V^gC 2 (IR 3 ) and thus A£ gC 2 (_B 3 ). 



Denoting by H V( the Hessian of Vg we obtain for each x G B 3 a point z G 0, x 
such that 



MKZ)(x)\< 



US(\*\) 



\Vt{x)-Vt 



+ 



tf<5(M)l 



\Wt(x) 



< 



u&\x\y 

( 'H v Jz)x,x) +^-\WV i {x)\<C\\D 2 V i \ 



\x 



\x\ 



Finally, for x&Si we have 



V(K0(tx) 



1 



(t)a ^(to)-V € (0)) + ^W € ( te ) 



XX. . . . t VVcCtx) 
x^^{Hv £ {rx)tx,tx) + — 



Ufa) 



t 



as t\0, uniformly in 16S1. We have shown that A£ G A, and since geX 
as seen above this implies that £g A. This completes the proof that Lo is 
onto and thus also the proof of Proposition [£T[ • 
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